A mechanical and simplified model for RC elements subjected to combined shear and axial tension

Very little is known about the shear behavior of elements, in particular those subjected to axial tension. The shear accompanied by tensile forces could cause premature failure of reinforced concrete, which is sudden with minimal warning. Therefore, understanding the shear behavior of reinforced concrete (RC) elements, including those subjected to axial tension, is an ultimate goal of the worldwide research community. In the current study, a new shear mechanical model for RC elements subjected to axial tension is developed, which makes physical sense and explains the behavior. The model is strain-based, inspired by the critical crack theory model (CSCT). In addition, the proposed model extended CSCT (ECSCT) quantifies the effect of axial tension forces on the shear strength in terms of reduction in the compression zone depth and increase in the longitudinal strain. Moreover, the nonlinear trend observed in the literature was implemented using nonlinear multi-variable regression. The ECSCT is validated and compared with available design methods with respect to an extensive database, including 180 elements tested under shear and tension from 18 different research investigations. The ECSCT provided an accurate and physically sound model yet safe to an acceptable extent. Last but not least, a simplified model for the purpose of design is proposed. The simplified model was chosen based on the mechanical model and calibrated using the extensive experimental database. The simplified model provided an accurate and simple model, yet safe to an acceptable extent.


Recent findings
For the last seven decades, several pioneering studies 26,[29][30][31][32][33] have been conducted to investigate the shear strength of elements under axial tensile forces; a brief recount of the most recent findings is as follows.
• The significant variables that affect the shear strength are element dimensions, reinforcement configuration, loading configuration, and boundary conditions; thus, the conclusions of different studies were inconsistent with each other 31 . For example, the device used to apply the axial tension could cause accidental restraint at the ends of the tested element 33 . • The angle of inclination of shear cracking is significantly affected by the axial tension, which makes it steeper 33 . On the other hand, the shear strength of beams with well-detailed longitudinal reinforcements is not affected by axial tension. This is due to the aggregate interlock mechanism. In addition, the compression longitudinal steel reinforcements decrease the effect of axial tension forces on shear, if any 29 . • A nonlinear relationship between the axial everyday tensile stresses and the shear strength was found by several researchers 32 . This contradicts the long-standing linear relation implemented by both the ACI 34 and the EC2 35 . • Pham 31 observed a decrease in the compression zone depth with the increase in the axial tensile forces, as shown in Fig. 2. • In general, design codes, including but not limited to the ACI and the EC2, are overly conservative, especially for cases of high tensile forces 26 . • An extensive experimental database of elements tested under shear and tension was gathered, combining the database complied by Deifalla 26 and Ehmann 30 . A total of 180 elements from 17 different research investigations. The data covered a wide range of all influential variables, as shown in Table 1 and Fig. 3 [30][31][32][33][36][37][38][39][40][41][42][43][44][45][46][47][48][49] . The effective parameters included the axial tension, the size, the shear-span to depth ratio, the concrete strength, flexure reinforcement ratio, and the width to depth ratio were gathered.

Model development
Introduction. This work is inspired by the Critical Shear Crack Theory (CSCT) 28 , which was first introduced in the 1990s and later implemented in the swiss design code, MC, and the new draft of the Eurocode 26,50 . The shear strength ( v u ) is calculated such that: while ω is the crack width, f ′ c is the cylinder compressive strength and d dg is the maximum nominal aggregate size.
(2) ω ∝ εd where d is the effective depth, and ε is the longitudinal strain, which is taken at 60% of the effective depth from extreme compression fibers 28 . In addition, the following assumptions were implemented, which are similar to the work by Deifalla 24,25 for slabs under combined punching shear and tension: (1) Plane cross-sections before deformation remain plain after deformation,m while maintaining small deformation.   www.nature.com/scientificreports/ distribution diagrams for a) flexure only, b) tensile forces only, and c) both, respectively; thus, the tensile strain at 60% of the effective depth from extreme compression fibers (ε) , which is calculated such that: www.nature.com/scientificreports/ where b is the element width, ρ is the flexure reinforcement ratio, N is the axial force ( positive is tension and negative is compression), M is the bending moment at the critical section for shear, E s is the steel reinforcements young's modulus (210,000 MPa), E c is the concrete young's modulus (10,000 3 f ′ c ), c is the compression zone depth, which is calculated such that: Figure 5 shows the variation of the compression zone depth (Eq. 4) versus the tensile stress for different flexure reinforcement ratios. Figure 5 was based on f Proposed failure criteria for combined shear and tension. Based on data shown in Table 1 and Fig. 3 as well as the failure criteria of the CSCT, the following form was proposed: Using nonlinear multi-variable regression, thus, the shear stress is calculated such that:     Table 3 shows the parameter table, including the probability of each parameter. It is worth noting that the power coefficient of the variable N/(bdρf y ) failed the hypotheses test, and it was found to be insignificant. Therefore, it is proposed that shear strength is calculated such that:

Models validation
For simplicity, two design codes were selected for comparison: the ACI and the EC2. However, it is worth noting that there is other design model that are more accurate with various levels of approximation, for example fib model code 50 . For simplicity, the model code was not selected as it requires a detailed calculation compared to the ACI and EC2. The strength calculated using the ECSCT and the simplified model were assessed against those calculated using the selected design codes with respect to the experimentally measured strength.
Several types of figures were implemented to compare the performance of the proposed models with the selected design codes, which were implemented in several investigations 12,51-54 . Firstly: a scatter plot between the measured and calculated strength was plotted for all models, which was assessed using the ideal 45-degree line and the inverse of the slope of the best-fitted line. While the strength in terms of stress is taken as the ratio between the shear force and the concrete cross area. Secondly: a histogram figure for the distribution of the ratio between measured and calculated strength (SR), which is assessed based on the distribution and being far from the ideal ratio of unity and lower coefficient of variation. The unity value for SR indicates the closeness of the calculated value to the measured one (i.e., the model accuracy). At the same time, the coefficient of variation of the SR distribution indicates the consistency of the model. In addition, the lower 95% of the SR indicates the safety of the model. It is the minimum SR value obtained using the model with a 95% confidence level. Therefore, the higher value of the lower 95% limit above the safety factor of design codes (approximately 0.85), the safer the model is. The confidence interval is calculated assuming a standard normal distribution. In addition, a significant level value of 0.05 represents the 95% confidence level. Thus, the lower 95% confidence limit is calculated using the following expression:   www.nature.com/scientificreports/ Thirdly: a scatter plot for the SR against various effective parameters is plotted, which is assessed using the inclination of the best-fitted line and the correlation coefficient (r). The closer the slope to zero is the more negligible effect of the variable on the accuracy and safety of the model. The r is the degree of association, which is such that: where, x and y are average values of variables x i and y i , respectively, and n is the number of tested specimens. The correlation coefficient is measured on a scale that varies from ±1 to 0. For example, ±1 , ±0.70 , ±0.50 , ±0.30 , and zero indicate exact, strong, moderate, weak, and no dependence, respectively. Thus, if the r value is less than ±0.30 shows that the model captured the effect of such parameter, while a coefficient ranged between ±0.50 and ±0.30 indicates a need for refinements in the modeling of this parameter. It is worth noting that the correlation coefficient is not directly related to the data scattering. Because the data scattering is dependent on the overall effect of the considered parameters in the specific model, while the correlation coefficient is indicative of the relation between the accuracy and a specific parameter. It is an indication not conclusive depending on the value as mentioned before.
Overall. Figure 6 shows the calculated shear strength versus the measured ones for the ECSCT, the simplified model, the ACI, and the EC2. While the strength is calculated in terms of stress taken as the ratio between the shear force and the concrete cross area. In addition, the line represents the actual performance and the linear fitted line for the model performance. Moreover, the inverse of the best-fitted line slope ( χ ) is indicated in the plots. The closer this value to unity is, the better accuracy and less divergence. The χ value for the ECSCT, the simplified model, the ACI, and the EC2 is 0.99, 0.88, 0.50, and 0.65. Thus, the strength calculated using the ECSCT and the simplified model is significantly less scattered than using the ACI and EC2.
Moreover, Fig. 6 shows the histogram of the SR values calculated using the ECSCT, the simplified model, the ACI, and the EC2. The SR is calculated using the simplified model, and the ECSCT is normally distributed around the ratio of unity compared to that using the ACI and the EC2. Last but not least, Table 4 shows the average, coefficient of variation, and lower 95% limit for SR calculated using the ECSCT, the simplified model, the ACI, and the EC2 for each study. It is clear that the performance of the ECSCT and the simplified model is more accurate and consistent than that of the ACI and EC2. However, it is safe with a lower 95% value of 0.96, higher than 0.85 targeted by most design codes. The simplified model and ECSCT SR values have a narrow range compared to that calculated using the ACI and the EC2, as shown in Table 4. ′ c , f y is the yield stress of the steel reinforcements. From Figs. 7, 8, 9 and 10, the safety of the strength calculated using the ECSCT and the simplified model is more consistent with the axial tension's effect than the ACI and EC2. The correlation coefficient (r) was calculated as 0.51-0.67, 0.54-0.73, 0.02-0.14, and 0.04-0.25 for the ACI, the EC2, the ECSCT, the simplified model, respectively. In addition, the slope of the best fit line for the SR calculated using the ACI, the EC2, is much higher than that calculated using the ECSCT, the simplified model. Thus, it is clear that SR calculated using the ECSCT and the simplified model are weakly correlated to the axial tension, while the ACI and EC2 are highly correlated . The variation of the tensile axial force does not affect the ECSCT and the simplified model compared to the ACI and the EC2 with respect to the experimental database.

Effect of shear span to depth ratio (a/d).
The SR is plotted against the specimen size in the shear span to depth ratio (a/d), as shown in Fig. 11. From Fig. 11, the safety of the strength calculated using the ECSCT and the simplified model is more consistent with the effect of the a/d compared with the ACI and EC2. However, the safety for the ECSCT is higher for non-slender elements with a/d value less than 3. This is because the original CSCT model was not developed for non-slender. In addition, the correlation coefficient (r) was calculated as 0.46, 0.49, 0.40, and 0.09 for the ACI, the EC2, the ECSCT, the simplified model, respectively. Thus, it is clear that SR calculated using the simplified model is weakly correlated to the a/d, while the ECSCT, the ACI, and EC2 are highly correlated to the a/d.
Moreover, the slope of the best fit line for the SR calculated using the ACI, the EC2, the ECSCT, and the simplified model is 0.46, 0.56, 0.10, and 0.02, respectively. The slope for the SR calculated using the ECSCT and simplified model is significantly lower than that calculated using the ACI and the EC2. The simplified model showed quite an improvement with respect to the effect of the arch mechanism in terms of the shear span to depth ratio (a/d).

Effect of specimen size (d).
Plots of the SR versus the specimen size in terms of effective depth (d) are shown in Fig. 12. From Fig. 12, the safety of the strength calculated using the ECSCT and the simplified model is more consistent with the effect of the d compared with the ACI and EC2. In addition, the correlation coef- Scientific Reports | (2022) 12:7863 | https://doi.org/10.1038/s41598-022-11577-y www.nature.com/scientificreports/ ficient (r) was calculated as 0.26, 0.32, 0.04, and 0.25 for the ACI, the EC2, the ECSCT, and the simplified model, respectively. Thus, it is clear that SR calculated using the simplified model, the ECSCT is weakly correlated to the d. Moreover, the slope of the best fit line for the SR calculated using the ACI, the EC2, the ECSCT, and the simplified model is 2300E-6, 2900E-6, 70E-6, and 1E-6, respectively. Thus, the slope for the SR calculated using the ECSCT and simplified model is significantly lower than that calculated using the ACI and the EC2. Therefore, the ECSCT model and the simplified model account for the effect of size in terms of d much better than the ACI and the EC2.
Effect of flexural reinforcement ratio ( ρ). The SR is plotted versus the size in terms of the flexure reinforcement ratio ( ρ ), as shown in Fig. 13. From Fig. 13, the safety of the strength calculated using the ECSCT and www.nature.com/scientificreports/ the simplified model is more consistent with the effect of the ρ compared with the ACI and EC2. In addition, the correlation coefficient (r) was calculated as 0.08, 0.00, 0.18, and 0.06 for the ACI, the EC2, the ECSCT, and the simplified model, respectively. Thus, it is clear that SR calculated using the ACI, the EC2, the ECSCT, and the simplified model is weakly correlated to the ρ . Moreover, the slope of the best fit line for the SR calculated using Effect of concrete compressive strength (fc′). The SR is plotted versus the size in terms of the concrete compressive strength (fc′), as shown in Fig. 14. From Fig. 14, the safety of the strength calculated using the ECSCT and the simplified model is more consistent with the fc′ effect compared with the ACI and EC2.  www.nature.com/scientificreports/ In addition, the correlation coefficient (r) was calculated as 0.08, 0.04, 0.04, and 0.01 for the ACI, the EC2, the ECSCT, and the simplified model, respectively. Thus, it is clear that SR calculated using the ACI, the EC2, the ECSCT, and the simplified model is weakly correlated to the fc′. Moreover, the slope of the best fit line for the SR calculated using the ACI, the EC2, the ECSCT, and the simplified model is 103E-4, 96E-4, 3E-4, and 11E-4, respectively. Thus, the slope for the SR calculated using the ECSCT and simplified model is significantly lower than that calculated using the ACI and the EC2. Therefore, it is clear that the ECSCT and the simplified model are more consistent and accurate with respect to the direct shear mechanism and the residual tensile stresses in terms of the concrete strength (fc′).

Effect of width to depth ratio ( b/d).
The SR is plotted versus the size in terms of the width to depth ratio ( b/d ), as shown in Fig. 15. From Fig. 15, the safety of the strength calculated using the ECSCT and the simplified model is more consistent with the effect of the b/d compared with the ACI and EC2. In addition, the correlation coefficient (r) was calculated as 0.29, 0.24, 0.29, and 0.22 for the ACI, the EC2, the ECSCT, and the simplified model, respectively. Thus, it is clear that SR calculated using the ACI, the EC2, the ECSCT, and the simplified model is moderately correlated to the b/d . Moreover, the slope of the best fit line for the SR calculated using the ACI, the EC2, the ECSCT, and the simplified model is 17E-2, 15E-2, 4E-2, and 2E-2, respectively. The slope for the SR calculated using the ACI and the EC2 is significantly higher than that calculated using the ECSCT and simplified model. Therefore, it is clear that the ECSCT and the simplified model are more consistent with respect to the aspect ratio.

Conclusions
A physically sound mechanical model capable of accurately reproducing the actual behavior of reinforced concrete members under combined shear and tension is proposed. In addition, the model is accurate and simple for design. The effect of axial tensile forces on the compression zone depth and longitudinal strain and, ultimately, on the shear strength is accounted for. The proposed model is based on the principles of mechanics and its applicability to reinforced concrete elements under shear combined with tension. In addition, a simplified model is proposed for the purpose of design. The two proposed models were found to be more accurate, consistent, and reasonably safe compared to selected design codes. Moreover, the effect of basic parameters on the safety of the proposed models and the selected design codes was assessed. For all basic variables, including (1) the axial tension; (2) the shear span to depth ratio; (3) the flexure reinforcement ratio; (4) the concrete compressive strength; and (5) the width to dept ratio, the following conclusions were reached.
-The correlation relation with the safety factor calculated using the American and European design code was very strong, while that for the proposed models was very weak. Thus, the proposed models captured the effect of all basic variables much better than the selected design models. -The slope of the best fit line for the safety factor calculated using the proposed models is very small compared to that using the selected design codes. Thus, the proposed models are more consistent with the basic variables than existing design codes.